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I. INTRODUCTION 



Although Einstein's gravitational field equations are in remarkable agreement with all 
solar-system and binary-pulsar tests they lead to a cosmo logical model which needs that 
the energy density of the Universe is strongly dominated by components, dark matter and 
dark energy, which have so far eluded direct detection. In particular, there is not enough 
baryonic matter to explain the observed properties of galactic dynamics using standard 
Einsteinian (or actually Newtonian) gravitational equations. The usual solution to this 
problem is to suppose that the vast majority of nonrelativistic matter in the universe consists 
of some weakly interacting particle we have not yet detected. Although there are several 
reasonable candidates for what this dark matter might be (see, e.g., [2|), it is worthwhile 
considering the alternative: It is possible that the field equations break down at galactic 
scales, i.e., that gravity is modified at distances relevant for dealing with galactic and inter- 
galactic dynamics. 

Along this line, Milgrom proposed a simple phenomenological law jsj] which leads to 
successful explanations of various observations, and to predictions 0] which turned out to 
be confirmed a posteriori. Milgrom's proposal, Modified Newtonian Dynamics (MOND), 
stipulates that a test particle at a distance r from a mass M will experience a gravitational 
acceleration given by the Newtonian expression = GM/r 2 as long as a at is (much) larger 
than a critical acceleration ao, while the same particle will undergo the MOND acceleration 
«mond = y/cLNdo = V GMclq / r when a at is smaller than a . (Although Milgrom's proposal 
can be viewed as a change in Newton's 3rd law (B| it is more often imagined as a modification 
of gravity, which is the view we shall take.) It turns out that the value 

a ^1.2xl0" 10 m.s" 2 (1) 

allows an excellent fit of galaxy rotation curves using reasonable mass-to-luminosity ratios 
Q, without the need for non-baryonic dark matter [!§]. MOND's ability to explain certain 
observed regularities of galactic structure contrasts favorably with dark matter, for which 
these regularities must either be accidental or else the result of some yet-to-be-discovered 
attractor solution in structure formation. Indeed, in the case of rotationally supported 
systems, MOND provides a simple explanation for (i) the Tully-Fisher relation [9j], which 
states that the observed limiting rotation velocity of galaxies, v^, scales as the fourth root of 
the baryonic mass of the galaxy (see [l(| for a recent dramatic confirmation of this relation); 
(ii) Milgrom's law, stating that the need for dark matter always seems to occur when the 
gravitational acceleration falls to about ao [HI; ( m ) Freeman's law, namely that the surface 
density never exceeds a /G [l2j; and (iv) Sancisi's law, i.e., that bumps in the rotation 
curves are correlated to the baryonic mass For pressure-supported systems, MOND 

also explains their typical size R = ■sj GMj a and predicts a stellar velocity dispersion 
a oc (GMao) 1 / 4 (l3l explaining thus the Faber- Jackson relation [l5j]. MOND also was 
able to predict properties of low surface brightness galaxies that were eventually confirmed by 
observation [lg, |l7| . Recently, Ref . [18] used a large catalog of widely separated binary star 
systems [l9[ as evidence for the breakdown of Kepler's third law at the MOND acceleration 
scale ao- 

On the other hand, the original MOND formulation does have some difficulties. In 
particular, a single acceleration scale ao fitted to galaxy rotation curves does not account for 
velocity dispersions in the cores of galactic clusters, which still need some amount of dark 



matter 20]. Similarly, X-ray and weak lensing data from the bullet cluster [2l| indicate 
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that dark matter exists at a different location from the gas. (Because weak lensing data is 
involved, one must make some assumption about how lensing occurs in MOND gravity, and 



this might affect the negative conclusion [2214241] .) A reasonable fraction of dark baryons 



and/or massive neutrinos is thus still required at cluster scales, even if the MOND scheme 
happens to be an actual law of Nature at galaxy scales. This would not be in contradiction 
with observation, since most of the baryons in the Universe have not yet been detected. Dark 
baryons are actually also required to explain the observed peaks in the CMB spectrum. 

The most serious problem of MOND is that it is not a complete theory, so that testing it 
often implies making guesses about its predictions for lensing or for cosmological evolution. 
This has not been for a lack of efforts, and a large number of theoretical constructions 
have been proposed over almost three decades to promote MOND to a consistent relativistic 
field theory. One major problem has always been simultaneously reproducing the Tully- 
Fisher relation and giving a sufficient amount of weak lensing. That problem was finally 
surmounted in 2004 by the tensor-vector-scalar (TeVeS) model constructed by Bekenstein 
(after years of work with Milgrom and Sanders) [25l-l30l|. in which the MOND force (implying 
the Tully-Fisher relation) is mediated by a scalar field, and where the presence of a unit 
timelike vector field helps in obtaining the right amount of light deflection from galaxies and 
clusters. The model has been shown to give better agreement with cosmological data than 



many believed any relativistic extension of MOND could do [31H37 . 

In its original formulation |28|, TeVeS suffers from several theoretical and experimen- 
tal difficulties 

notably a serious instability |4ll 42]. The latest version of TeVeS 



(43I . 44]|, inspired by the Einstein- Aether framework (45 -50], seems to avoid this instability 
and predicts post-Newtonian parameters consistent with solar-system tests. However, it 
still needs an unnaturally fine-tuned function of the scalar kinetic term in its action to be 
also consistent with binary-pulsar tests j38[. (The extended Vainshtein mechanism recently 
proposed in [5l[ is a way to avoid this difficulty.) Out of the many alternative models which 
have been proposed in the literature (for example, see (H3|), the recent bi-metric theory [53[ 
is a particularly promising and elegant one, although its detailed properties (notably its 
stability) remain to be fully understood. 

As promising as we consider TeVeS to be, its dependence on other fields to carry part 
of the gravitational force is somewhat counter to the spirit of relativity. In the present 
paper, we re-examine pure-metric formulations of MOND along the general lines previously 
considered in [54] Our aim is not to here produce the ultimate theory but rather just to 
show what form any pure metric generalization of MOND must take in order to combine 
two key features for a static, spherically symmetric and pressureless source which contains 
no dark matter: 

• reproduce the MOND force law in the ultra-weak field regime of accelerations compa- 
rable to ao; and 

• produce enough weak lensing to be compatible with observations. 

We first derive the form the MOND corrections to the Lagrangian must take in order to 
combine these properties when specialized to a static and spherically symmetric geometry. 

1 By "pure-metric" , we mean that the full gravitational interaction is described by the dynamics of a single 
metric tensor g^ u , without introducing explicit extra fields like scalars or vectors, although such degrees 
of freedom may actually be hidden in some excitations of g^ v . For instance, we would call pure- metric the 
class of f(R) models, although it is well known they are equivalent to specific scalar-tensor theories. Our 
phrase pure- metric should also be distinguished from what is called a "metric theory" in [l| , meaning there 
that matter is minimally coupled to a single metric tensor g^. What we call pure- metric is a subclass of 
such metric theories, but we also impose that the^inetic term of gravity itself is a functional of only g^ v . 



Then we demonstrate that no local curvature scalar has this form. However, nonlocal 
scalars do exist which take the correct form, and we exhibit some. As anticipated in |54j |. 
the Lagrangian inevitably becomes cubic in the weak fields, raising concerns about stability 
which we discuss briefly in the conclusion. 

This paper is organized as follows. In Sec. [Til we discuss some basic phenomenological 
properties used to deal with static, spherically symmetric systems, and we show how a pure 
metric action reproducing the MOND phenomenology can be devised for such systems. In 
Sec. IHIl we demonstrate that such an action cannot be local, i.e., it cannot be a function 
of only the metric and a finite number of its derivatives. Section [IV] introduces the main 
ingredients needed to construct a suitable nonlocal action for gravity. In Sec. [V]we exhibit a 
nonlocal model having the properties discovered in Sec. [UJ i-e., which reproduces the MOND 
dynamics at large distances, including enough weak lensing, while tending towards general 
relativity at small distances. Our conclusions are given in Sec. [VI] 

II. PHENOMENOLOGY 

The point of this section is to derive the form that the MOND modification to the 
gravitational Lagrangian must take when specialized to the ultra-weak field regime of a 
static, spherically symmetric geometry, 



We begin by reviewing how the equations of general relativity work for a source which would, 
of course, need to consist mostly of dark matter. In the ultra- weak field limit, these equations 
imply relations for the two linearized potentials, a(r) = A(r) — 1 and b(r) = B(r) — 1. One of 
these relations determines how the potentials depend upon the source and the other fixes how 
they depend upon each other. Our metric interpolation of MOND consists of changing how 
the potentials depend upon the source but not much how they depend upon one another. As 
the section closes we consider the form the MOND correction to the gravitational Lagrangian 
must take in order to substitute our MOND equations for those of general relativity. 
We assume a perfect fluid source, 



where p(r) and P(r) are respectively the energy density and pressure. Only two of the ten 
field equations are independent in this geometry; the rest are either trivial or implied by 
conservation. Defining G as Newton's constant, the tt and rr Einstein equations equations 
are, 



ds 2 = -B(r)c 2 dt 2 + A(r)dr 2 + r 2 dQ 2 . 



(2) 



T; = diag(-p,P,P,P), 



(3) 



G, 



tt 



A' 




(4) 



G 



B 



rr 



rA 2 

B' 



(5) 



.4 



rAB 



Equation ([4]) can be integrated to give us the rr component, 



A(r) 



1 



2GM(r)i 



(6) 



c 2 r 
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where the enclosed mass is 



47T 



M(r) = ^ I dr'r' 2 p{r'). (7) 
c Jo 

The second equation (j3j) could also be integrated but we shall not need to do this. 

Now consider the regime of zero pressure and very weak potentials, for which the linearized 
potentials take the form, 

. . 2GM(r) 

a(r) « — -^-L m rb'(r). (8) 

These relations can be expressed in many ways but a convenient form, for our purposes, is 
as one equation for how the potentials depend upon the source, 

llt , 2GM(r) 
rb'(r) « — ^ , (9) 

and another equation for how the two potentials depend upon each other, 

a(r)^rb'(r). (10) 

The first equation (Q is what tells us that explaining cosmic motions requires dark matter, 
whereas the second equation (flOl) tells us that the amount of weak lensing is consistent with 
the data, assuming cosmic motions are explained. 

For circular geodesic motion at fixed radius r with angular velocity 0, one can show 

2r 2 6 2 2v 2 

rB'ir) = rb'ir) = = Ar . (11) 

c 2 c 2 

We emphasize that relation (flT|) depends only upon the geometry ([2]) and minimal coupling 
to matter (that we will always assume within the present paper), without any assumption 
about the gravitational field equations which produce it. The Tully-Fisher relation implies 
that the rotational speed v (r) = r<p tends to a constant which goes as the fourth root of the 
source luminosity. MOND imposes the Tully-Fisher relation by changing equation to [8J , 

c 2 

where the right arrow indicates that the relation applies in the ultra-weak field regime of 
low accelerations. 

Relation (fl2|) contains the physics we want, but it is not yet in the form of a modification 
to just the left hand side of the gravitational field equations (jUl5]). To reach that form 
we need to isolate the local energy density p(r) by first squaring, then differentiating and 
shifting some factors from right to left, 

The other MOND equation can be written in a variety of ways because the right hand side 
vanishes for the relevant case of zero pressure. The weak lensing data is also not good enough 
to justify insisting upon precisely (ITU1) . so we would be happy with a{r) = krb'(r) for any 
positive, order one constant k. This suggests the second MOND equation should take the 
form 

c 2 / \ 2 

-(krb'-a) =0. (14) 

a r d V / 
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Note that it would not change the MOND phenomenology were we to multiply f lT4|) by a 
constant; we could also add a constant times it to ( TT3|) . 

Relations ( TT3lfT4|) are the modified gravity equations we wish to attain in the ultra-weak 
field regime for a static, spherically symmetric and pressureless source. We now seek an 
ultra-weak field expansion of a Lagrangian £mond which cancels that of general relativity 
£eh and substitutes cubic terms whose variation gives (I13til4|) . Although one generally loses 
field equations by specializing the metric before variation, we shall recover the correct g u 



and g rr equations |55|, |56j], in the ultra- weak field regime of course. The equations lost by 
specializing first are those associated with conservation. 

After some judicious partial integrations, the Einstein-Hilbert Lagrangian takes the form 



£eh = — „ R\Z^g — > Surface term) + —{ -rab' + — + 0(h 3 ) \, 
1d7tG V / lb7rG [ 2 J 



(15) 



where h stands for a and b. In the ultra-weak field regime, the MOND Lagrangian we 
seek should have a quadratic term that cancels the quadratic part of the Einstein-Hilbert 
Lagrangian, plus a cubic term which enforces our interpolation ( JT31TH1) of the MOND physics. 
The most general Lagrangian of this form is 



^MOND _ 



16ttG 



a 



2 



2 



c 



era 



3 



rab' + 0{h i ) + + Pa 2 b' + 1 rab' 2 + 5r 2 b' 3 + 0{h i ) }, (16) 



a l r 



where a, /3, 7 and 5 are dimensionless constants whose properties we shall constrain using 
the phenomenology of the ultra-weak field regime. A minor point which deserves comment 
is that the 0(h 3 ) corrections to the first square-bracketed expression in ffTBT) differ from the 
cubic MOND terms in the second square-bracketed expression by a factor of aor/c 2 , which 
would only become of order one on horizon scales and is utterly negligible on galaxy scales. 

The gravity Lagrangian is £g rav = £eh + £mond, anc ^ we w i sn to compute the variation 
of the associated action when specialized to a static, spherically symmetric geometry. Recall 
that the full Einstein equations for arbitrary geometry are obtained by varying the Einstein- 
Hilbert action as 

16nG 5S EH 87rG 

G^{x) = — -Tp,{x). (17) 



c 4 y/=g5g^(x) * w c 4 

For a static, spherically symmetric geometry \J—g = r 2 ^J (l + a)(l + 6). We want g tt = 
— (1 + a) and g rr = 1 + 6, and we neglect higher powers of the weak fields, so the relevant 
variations for us (assuming zero pressure) are 



16nG 5S grSkV 


c 2 


c 4 r 2 £5( r ) 


aor 2 


16nG 5S graiV 


c 2 


c 4 r 2 5a(r) 


aor 2 



I3{a 2 )' + 2 7 (ra6 / )' + 35(r 2 6' 2 )'| = , (18) 

^! + 2 /3a6' + 7 r6 ,2 l = 0. (19) 



Demanding that equation ffl9l should have the unique solution a = krb' implies 

a = - A an d 7 = -/3k. (20) 

OrC 

Substituting ff20l into equation fflBl and demanding that it give f|T3l) implies 

6 = -l + l^ 2 - (21) 
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Hence the MOND Lagrangian we seek has the following expansion in the ultra-weak field 



regime: 
£mond 



c 4 r 2 
16nG 



\ab' a 2 ^.,o 
V-2? + ° {h 



C 2 ' 
+ — 

a 



where the constant (3 must be nonzero but is otherwise arbitrary. 



III. LOCAL TOOLS FOR MODEL BUILDING 



(22) 



In the previous section we considered static, spherically symmetric geometries in the 
ultra-weak field limit for which MOND ought to apply. Our result is that the form ( 122]) 
for the MOND addition to the Einstein-Hilbert Lagrangian, allows us to reproduce the 
MOND force and sufficient lensing without dark matter. The burden of this section is that 
no local, invariant Lagrangian can have that form. Of course minus the Einstein-Hilbert 
Lagrangian reproduces the quadratic parts of (1221) . so it is the cubic terms which comprise 
the fundamental obstacle. Our analysis is a straightforward proof by exhaustion: we examine 
all local curvature scalars for static, spherically symmetric geometries in the ultra-weak field 
regime. 

It is useful to expand the various curvatures in powers of the graviton field h^ u = g^—f}^. 
We can greatly simplify the analysis by working in Cartesian coordinates, which makes the 
affine connection vanish in the absence of curvature. The 3 + 1 decomposition of the graviton 
field is 



*oo 



—b(r) 



h 



0, 







and 



hi 



a(r)Pr ? 



where r* = x % jr. It is also useful to introduce the projector 



7T- 



(23) 
(24) 



The nonzero components of the affine connection are 
b' 

r' + (){lr) , I '„, ■ 



r° 



oi- if + 0{h 



2\ -pi 
J 1 00 



^P + 0(/i 2 ) , ^^ fc = ^Pf 5 f* + ^p7^^ + 0(/^/^ , ). (25) 



In denoting higher order corrections, we make no distinction between derivatives and inverse 
powers of r. So the term O(hh') includes terms of the form hh! and h 2 jr. 
Our convention for the Riemann tensor is 



= f) yf — f) V p -4- V p V a —V p V a 



Its nonzero components for a static, spherically symmetric geometry are 



h" h' 



-r*'f J + —tt ij + 0(ti 
2 2r v 



2r 



a 

+ ^7 



—ik„ji „il„jk 
71 7T — 71 7T J 



0{h' 2 ). 



(26) 

(27) 
(28) 



We define the Ricci tensor as R 



R" 



Rqo 
R 



fipu ' 
2 X 



'J 



b" b' 

T + - + o(h' 

2 r 

— — H — Pr J + 
2 r J 



and its nonzero components in our geometry are 

(29) 



V a 1 a 
2r 2r r 2 . 



7T ij + 0{h 



>2n 



(30) 



7 



The Ricci scalar is R 



y 1 L fiu -i 
R = 



and it works out to be 
-b" 



2b' 2a! 2a ^., /2 , 

— + — + ^ + 0{h' 2 ). (31) 

ry ry ry £ 

From the preceding analysis we note that every nonzero component of the curvature 
involves two derivatives (or inverse powers of r) acting on one or more weak field, 



Curvature ~ h" + 0{h' 2 ). 



(32) 



No matter how the indices are contracted, N factors of the curvature must therefore have 
the form 

(Curvature) " ~ (ti'f + O ((/OW* -1 ) • ( 33 ) 
we seek involves powers of just one derivative acting on a single 

„2 



The MOND correction 
weak field, 



C 



MOND 



c 4 r 2 
16ttG 



{h'f + -{h'f + 0{h A ) 
a 



(34) 



The Ricci scalar ( 13T|) gives the quadratic terms because its linear part is a total derivative. 
However, the cubic terms of (I34p not only have too few derivatives per weak field, they also 
contribute an odd total number of derivatives. The latter problem is much worse than the 
former because the leading weak field term in a curvature scalar might drop out — as it does 
for R — but nothing can change the total number of the derivatives it contains. Including 
differentiated curvatures increases the number of derivatives per weak field, and can in any 
case only add an even number of derivatives once all the indices are contracted to form a 
scalar. 

That completes the main argument of this section but it is worth giving the nonzero 
components of the Einstein and Weyl tensors for future reference: 



Goo 

GiA 



- + — + Q{h 



a 



OiOj 



b' 

r 

1 

~Y2 L 
1 



j^lyO _|_ 



2 



2r 



a 
2r 



ir ij + 0(h' 2 



,„ V a' 2a 
b" - - + 



<) ,J 3r'rM +0(h 



,/2-n 



ijki 



6 



b" 



b' 

r 



a 
r 



2a 



\ Lhj ry^ry 1 ^ 

2 



2 



2 



l2s 



(35) 
(36) 
(37) 

(38) 



Note that all components of the Weyl tensor are proportional to the same linear combination 
of the weak fields, so that any scalar formed from C pcTUV will access this combination 

1 

3 



^Jp(T(J,U £1 



popv 



-b" + -- 



b' a! 2an 2 



+ 



+ 0{h"h' 2 ) 



It is also worth noting some of the other scalars we can get: 



R 1 



,.. 2b' 2a! 
-b" + — + — 



R pafll/ Rpa-fiv — ^R^Ruu + R 



2 



2a- 


2 






4 / 




— ab 







+ 0(h"h 



[ab" + a'6') + 0{h"h' 2 " 



(39) 

(40) 
(41) 



8 



IV. NONLOCAL TOOLS FOR MODEL BUILDING 



The new features that nonlocality brings to model building are that inverse differential 
operators reduce the number of derivatives, and that the gradient of the invariant volume of 
the past light-cone allows us to define a timelike 4- vector with which we can select particular 
components of the curvature. The first feature is necessary because, as discussed in the 
previous section, curvature scalars involve powers of two derivatives (or factors of 1/r) 
acting on a weak field, whereas the MOND correction (1221) we seek to realize as a scalar 
involves powers of only a single derivative of a weak field. The second property is needed to 
get the right weak fields. 

A philosophical digression is necessary at this point. We do not maintain that physics is 
nonlocal at the fundamental level; we believe rather that nonlocality enters through quantum 
corrections to the effective field equations from loops of massless gravitons. These induce no 
macroscopic nonlocality in flat space background because their interactions are suppressed 
by derivatives, however, the situation is quite different when a cosmological constant is 
present. It has been argued that self- interactions between elements of the vast ensemble of 
infrared gravitons produced during primordial inflation show secular growth which eventually 
becomes nonperturbatively strong 57]. Nonlocal effective field equations for cosmology have 
been studied [58| as a way of abstracting these effects to the nonperturbative regime. Our 
work here will apply the very same nonlocal tools to build a model of structure formation. 
Although we work on a purely phenomenological level, it might be possible to derive a 
successful model from first principles using the Schwinger-Keldysh formalism [59]. Hence 
our nonlocal constructions will always be viewed as proceeding from causal evolution, based 
on the notion that the universe was released in a prepared state at some finite time. This 
last point is the key to being able to define a timelike 4- vector field and it must be accepted, 
even if one chooses to disregard our motivations and treat the models we propose on a purely 
phenomenological level. 



A. The inverse scalar d'Alembertian 



The scalar dAlembertian is familiar to students of general relativity, 

1 



□ 



--d„ 



'-9 



(42) 



We define the function F(x) obtained by acting □ 1 on any function f(x) (that is, F(x) 
D^ 1 /) by the solution of the differential equation 



DF(x) = f(x) 



(43) 



subject to retarded boundary conditions. Specializing to the case of a static, spherically 
symmetric geometry, and a source function f(r) which falls off at infinity, we obtain an 
equation which can be solved by integration, 



1 



d 



r 2 ^A{r)B{r) dr 



2 j Bjjr) dF(r) 
A{r) dr 



f(r) 



2 j B{r) dF(r) 
A{r) dr 



dr' r' 2 ^A{r')B{r') f(r') 



(44) 
(45) 
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F(r) 



^\l§0) f r dr"r"WA(r»)B(r»)f(r"). (46) 



In the weak field limit for a source f(r) which is already first order we can set A — B — 1. 
We can also change the order of integration and perform the r' integration to get 

i r°°dr' r' 

y=-j r dr " r " 2 f( r ") + °^ 2 ) ( 47 ) 

= - / dr'r' 2 f{r') X - - / dr'r' 2 f(r') x - f + 0(h 2 ) (48) 

JO r Jr r ' 

1 /> x 'p™ + 0(/i 2 ). (49) 



6*( ct — llx— x'l 



47T y 1 1 a? — af' 

Of course this is the usual Coulomb Green's function. We can make similar contact with the 
Lienard-Wiechert potential if we regard the system as released at some early time labeled 
t = 0, 

1 r r°° r s[c(t— t')— \\x— x'W ) 

□ / = -^ I dt'jd'x'^ ^-^1 J -xf(\\x'\\) + 0(h 2 ) (50) 

= ~ Jd 3 x' V p^jj x /(H* ||) + 0(/, 2 ). (51) 

The theta function in expression (|5ip is usually irrelevant for functions f(r) which fall 
off rapidly, and for late times t. However, it plays an important role when the function f(r) 
happens to be constant, 

l/o = -f j C 'dr'r> + 0{h 2 ) = -\h{ct) 2 + 0{h 2 ). (52) 

It is well to remember that even our static, spherically symmetric systems are embedded 
in a larger cosmological background which had a beginning and is even now slightly time 
dependent. 



B. A timelike 4- vector field 

The preceding considerations are especially important for our second nonlocal building 
block: the invariant volume of the past light-cone. Suppose S is the Cauchy surface on which 
the initial state was released and let M. stand for the spacetime manifold comprising S and 
its future. For a general metric g^ u we define the invariant volume of the past light-cone 
from the spacetime point as 

V[g)(x) = J^x'^gJ^)e^-a[g}{x,x'))e(^[g}{x,x')). (53) 

Here a[g](x,x') is the geodesic length function introduced by DeWitt and Brehme 60j|. In 
expression ( |53l) . we note 9{j r [g](x,x')), a functional which is the invariant generalization 
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of 6(x° — x'°) needed to restrict the integration over x" 1 to the past of x^. [T stands 
for "forward" in this notation.) This functional is defined as one when the extension of 
the geodesic between x^ and x'^ eventually intersects the initial value surface S, and zero 
otherwise. 

The volume of the past light-cone is of great interest to us because it is guaranteed to 
grow when the point x M evolves in whatever is the timelike direction of the metric q, iv . Hence 



its gradient must be timelike and can be used to define a timelike vector field [58 



g^{x)d v V[g]{x) 



u ft \a)(x) = 

yf-g*{x)d a V\g]{x)d p V\g](x) 

For the static, spherically symmetric geometry we have been considering it reduces to 



(54) 



u»[g](x) 



rS ft 



Br) 



(55) 



It can therefore be used to pick out the timelike components of a tensor, just like the 
fundamental vector field of TeVeS [28| . 

For our purposes it is better to exploit the close relation which exists between the volume 
of the past light-cone and the functional inverse of the Paneitz operator, 



D P [g) = D 2 + 2D, 



RT - -g^ u R 
3 y 



D„ 



(56) 



This 4th order differential operator appears in conformal anomalies [6jJ. The relation be- 
tween it and the volume of the past light-cone is that 871 /Dp acting on one agrees with V 
for arbitrary homogeneous and isotropic spacetimes [621 ]. 



57T 



Dp [FRW] 



1 = V[FRW]. 



(57) 



62| but that is probably irrel- 



Perturbations away from this background do not quite agree 
evant for any use we might make of V[g](:r). 

The great advantage to defining our timelike vector field using the inverse of Dp is that 



we can avail ourselves of a simple partial integration trick [63|, [64J for deriving causal and 
conserved field equations. To understand the trick, consider varying the product of a local 
functional of the metric F[g] times some inverse differential operator T>~ x — either or 
acting on another local functional G[g], 



Dp 1 



6 



5g» u (x) 



V 5gt* u V 



V5g 



(58) 



The second and third terms on the right of expression (1581) would make acausal contributions 
to the field equations which involve fields to the future of x* 1 . 

To see the acausality of expression ( 158]) more clearly, let us expand out the final term, 
with all the implied integrations and coordinate dependence made explicit. In order to fix 
notation we express the term being varied in (158|) as 



V 



d 4 x' F(x') d 4 x" g Tet (x'; x")G(x"). 



(59) 
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Here Qret(x'; x") is the retarded Green's function associated with the differential operator V, 
and "retarded" means that it vanishes for x"° > x'°. In this same language the final term 
on the right of (1551) would be 

Fi-^- = k'F(x') f^ x "g rct ( x ';x")^p- . (60) 
V8g» v J y J J y J 5g^(x) v ; 

Saying G[g\(x") is local means it depends only on the metric and some finite number of its 
derivatives at x' m . Hence its variation with respect to g tJ,y (x) is proportional to at most a 
finite number of derivatives of S 4 (x"—x). Of course this means we can perform the integration 
over x'" 1 to get at most some derivatives acting on Q ret (x';x) = Gad v (x;x'). The remaining 
integration over x" 1 involves fields to the future of x^. 

This sort of acausality is inevitable for any nonlocal action based on a single field. The 
Schwinger-Keldysh effective field equations avoid it by the same physical field being repre- 
sented in a complicated way with two dummy fields. One first varies with respect to one 
of the dummy fields and then sets the two dummy fields equal, after which cancellations 
between various contributions result in there being no dependence upon dynamical variables 
to the future of x^. We shall circumvent this complication by having recourse to the simple 
trick of "partially integrating" the acausal terms of (155|) so that their nonlocality is restricted 



to the past of x M [63|, [64 



F L™L G ^J™L G )L F , (61) 

VSg^V \8g^V JV v ' 

„ 1 5G 5G 1 

F. (62) 



V bg^ v 5g» u V 

The result is manifestly causal. It is also conserved (if we include the variation of the mea- 
sure factor) because we have just substituted, in the field equations, the causal retarded 
Green's function everywhere an acausal advanced Green's function appeared. Conservation 
requires only the differential equation, which both the advanced and retarded solutions obey. 
Of course this is just a trick; a true derivation from fundamental theory would require use 
of the Schwinger-Keldysh formalism [59]. However, the object of our study is the effec- 
tive field equations, and they are perfectly valid as long as we consider them on a purely 
phenomenological level. 



V. AN EXPLICIT MODEL 

There are many ways to define a suitable relativistic generalization of £mond- A partic- 
ularly elegant construction is based on two nonlocal building blocks, 



X[g](x) = g 



d^(R aP u a u^-^R 



d, J h(Rp (T u p u a -\R 



□ 



Y[g](x)=g^ d,-[2R aP u a u? 



O v -[-lR llo u' 



(63) 
(64) 



Although these scalars are deeply nonlocal, even when specialized to static and spherically 
symmetric geometries, they give local, and very simple results, to lowest order in the weak 
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field expansion. To derive these limits recall first the weak field, static and spherically 
symmetric results for i? o and R from expressions (129|) and (13T1) . respectively, 

Roo^^(r 2 b')' + O(h 2 ), (65) 

R — > (-r 2 6' + 2ray + 0(/i 2 ) . (66) 

The arrow indicates specialization to static, spherically symmetric geometries in the weak 
field limit. Note also that our statement that the residues are u O(h 2 )" refers only to their 
dependence upon the weak fields, without regard to derivatives or powers of r. 

The specialization of the 4- vector u^[g](x) to a static, spherically symmetric geometry is 
given by expression (155]) . Hence it is just u^[g\(x) — > 5q + 0(h) to the order we require, 
and we can write, 

R a pu a u fi --R^ ^{r 2 b'~raj + 0(h 2 ) , (67) 
«„/9 , jL/WY , n(h 2 



2R a pU a u p — >• — ( r 2 b' ) + 0(h 2 ) . (68) 



These terms are both first order in the weak fields so one can use expression (1481) to implement 
the action of D _1 on them, 

l[R apU ^- l -R) — ► - fdr'{b'(r')-^l) + 0(h 2 ) , (69) 

°,A _^ _ J dr' b\r') + 0(h 2 ) . (70) 



^{2R aP u a u 



Of course the only derivatives that matter are with respect to the radial coordinate r 

1 /_ „ .a 1 



d,-(R aP u^--R) -^%(y-*) +0(h 2 ) , (71) 
(2R a(3 u a u^ — ► ^ + 0(/i 2 ) . (72) 



The residue terms in these expressions are still nonlocal. However, substituting (ITTj) and 
f!72|) into expressions fl63l) and (164|) gives a local result to leading order, 

X[g]{x) ^ (b' ~)\ 0{h*) , (73) 
Y[g]{x)— >b' 2 + 0(h*). (74) 
Choosing k — 1 in f )22p . we find that acceptable MOND equations would result from 

£ ™ - t4{ K- A ' +y ) + ^(w 1 -^ 1 ) + ■ ■ W ■ < 75 > 



Relation ( 1751) gives just the first two terms in the ultra-weak field expansion of the the- 
ory. That is all we can infer from the deep MONDian regime. There are many ways of 
extending the expansion to all orders to define the full theory. The chief requirement on 
any such extension is that it be suitably suppressed in comparison with general relativity 
for Newtonian accelerations much larger than qq. Newtonian gravity seems to be valid in 
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the solar system out to at least 80 Astronomical Units (the furthest of the Pioneer Probes), 



at which point <7jv/oo ~ 10 [65(] . In this regime we can take the Newtonian acceleration to 
be gN ~ c 2 b'{r) ~ c 2 a(r)/r, so we need 



C 



MOND 



< 



c 



GR 



16nG \a ) 



(76) 



for g N /a > 10 4 . 

It is best to study the weak-field regime using dimensionless variables 



c 2 


&'- 


a 


3ao 


r 


c 2 








-C 


3ao 







0{h 2 



0{K 2 



(77) 
(78) 



Whereas the variable y[g] is of order one or smaller in the ultra- weak regime, it is of order 
Sn/O'O ~ 10 4 in the solar system. However, the variable x[g] vanishes, to lowest order, in 
both regimes, so whatever function interpolates between the two regimes must involve y[g]. 
With the variables ( J771) -(I78l) . the ultra- weak field expansion of the MOND Lagrangian 
f5l) takes the form 



C 



MOND 



9a 2 
32ttG 



-x 2 + y 2 + x 3 - y 3 + 



(79) 



We are therefore seeking an extension of the bracketed term in (IT9]) which is suppressed, 
relative to y 2 , for large y and x ~ 0, and whose corrections to y 3 are numerically small for 
y 1- Of course many functions of x and y have this property. However, we also need to 
pass tests of post-Newtonian gravity in the solar system and in binary pulsars, therefore the 
suppression of ( 1791) should be very efficient at small distances. An extra constraint on any 
possible extension of f )79p is that its variation with respect to x (i.e., to the radial component 
of the metric, a) should allow the looked-for solution x = (i.e., a = rb'). We just quote 
here two examples of such extensions having the required properties, and their associated 
behaviors for large y and x ~ : 



{y — x) x (^x + xy + y j < 



e » — x e 



y 2 e-y, 
y 2 e~ y2 - y . 



(80) 
(81) 



It is easy to check that the predicted deviations from general relativity are exponentially 
small with respect to the tightest solar-system constraints, but that the MOND behavior 
f!79|) is predicted at large distances. 

The MOND Lagrangian fl22|) was constructed in Sec. [Til in order to cancel the general 
relativistic predictions at large distances while imposing the precise physics we wished to 
reproduce. In particular, we saw that it was possible to predict any amount of weak lensing 
by changing the numerical value of the coefficient k. In the present section, we chose k = 1 
to recover the same weak lensing as predicted by general relativity in presence of a dark 
matter halo. In such a case, it is not necessary to cancel the x 2 term coming from the 
Einstein-Hilbert action and to add a cubic x 3 as in (!79|) above. Indeed, the original x 2 term 



is enough to force x 
instance 



0, and we may thus consider a Lagrangian depending only on y, for 



C 



MOND 



32ttG 



(82) 
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Added to the Einstein-Hilbert term, this suffices to reproduce the MOND dynamics and 
enough weak lensing at large distances, while predicting fully negligible deviations from 
general relativity at small distances. 

None of these Lagrangians (I8UH82P is analytic in a , but they all possess the key property of 
vanishing when ao goes to zero from above. To see this, note that they vanish for x = = y, 
irrespective of ao- Note also that neither x nor y can be negative, so if y = k/ao for some 
positive constant k, then the limiting form, for small ao, of flHOj) and fl82|) vanishes like e _fc//a °, 
while the limiting form of (ISTj) vanishes like e~ k / a °. 

Our final comment on explicit models concerns the "external field effect" in which MON- 
Dian behavior of one system can be severely affected by another [66]. This property is 
deeply embedded in the nonlocal constructions of our scalars X[g](x) and As one 

can see from their definitions (I63H64I) . these scalars involve the nonlocal operator D _1 acting 
on curvature scalars which are themselves contracted into the normalized gradient u^[g\(x) 
of the invariant volume of the past light-cone. In the static, spherically symmetric limit we 
have studied, X[g](x) and y[g](a;) depend only on the central gravitating source. However, 
they can be quite different, even in the static limit, when other sources are present. It is 
highly significant that they also depend upon past history. This holds out the possibility 
for reconciling problems in describing recently disturbed systems such as the Bullet Cluster 



|22l . 123] . Of course we cannot, at this stage, claim that our model incorporates the ex- 
ternal field effect in a desirable way; what actually happens beyond the static, spherically 
symmetric limit is a matter for future study. 



VI. DISCUSSION 

We have considered the problem of devising a pure metric interpolation of MOND, with 
neither dark matter nor additional fields, for static, spherically symmetric systems. In 
the deep MOND regime of small accelerations, gravity is described by two weak fields, 
6(r) = —gtt — 1 and air) = g rr — 1. In this regime the MOND force law is given by equation 
( fl2|) . and the requirement that there be enough weak lensing is roughly a(r) = krb'(r) for 
some positive constant k of order one. Our first result is that the ultra-weak field limiting 
forms of the g u and g rr equations are ( !T3irH|) . subject only to the ambiguity of multiplying 
( fl4|) by a constant or adding such a term to (Fl3|) . Our second result is that reaching this 
form requires the full gravitational Lagrangian £ grav to possess a MOND correction to the 
Einstein-Hilbert term, £ gr av = £eh + ^mond, where the ultra- weak field expansion of this 
correction takes the form f)22p . 

We then turned to how the MOND Lagrangian £mond depends upon a general metric. 
Our third result is that no local curvature scalar can reproduce the ultra-weak field form 
(|2"2"|) . The reason is that curvature scalars involve powers of two derivatives of a weak field 
whereas the MOND correction fl2"2l involves powers of only a single derivative of the weak 
fields. 

Nonlocal models have the great advantage that they allow one to effectively remove 
derivatives. Our fourth result is that it is possible to construct invariant nonlocal models 
which degenerate to (1221) . for static and spherically symmetric geometries in the ultra- weak 
field limit. In fact there seem to be many ways to do this, some of which are laid out in 
section |Vl So it would be fair, at this stage, to say we are developing a class of models rather 
than a unique model. 

As explained in section IIVt our constructions involve two nonlocal building blocks: the 
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inverse scalar d'Alembertian f H2|) and the timelike vector field u^[g\(x) formed from nor- 
malizing the gradient of either the volume of the past light-cone (|53|) or the closely related 
inverse of the Paneitz operator (157|) . Unlike TeVeS, the timelike vector field of our class of 
models is not an independent variable but rather a nonlocal functional of the metric itself. 
In our view this nonlocality is not fundamental but should be viewed rather as the result of 
quantum corrections (perhaps from the epoch of primordial inflation) to the effective field 
equations. So one should always bear in mind that our class of models involves the universe 
being released in some prepared initial state at a finite time. A derivation from fundamen- 



tal theory would be in the context of the Schwinger-Keldysh formalism [59j. In the purely 



phenomenological context of our current work, we employ the partial integration trick ( 162]) 



introduced in [63|, |64| to derive causal and conserved field equations. 



Although the timelike vector field u^[g\(x) will certainly introduce preferred frame effects, 
we believe these should only be significant in the ultra-weak field limit for which the MOND 
corrections become important. Even in this regime they should be suppressed by the square 
of a peculiar velocity divided by the speed of light. Typical peculiar velocities are several 
hundreds of kilometers per second, so the suppression factor should be about 10~ 6 , which 
is not likely to be observable. However, it may be very significant that our model depends 
upon events in the past light-cone. One consequence of this dependence is that recently 
disturbed systems such as the Bullet Cluster may be far from the static MOND limit. 

With any of the full metric interpolations described in section |VJ it would be possible 
to study the important issues of cosmological evolution and stability. As anticipated in 
our gravitational equations (I13til4|) are quadratic in the ultra-weak field regime, which 



means the gravitational Lagrangian is cubic. That poses an obvious potential problem for 
stability, although our fears on this score might be avoided by the absolute values needed 
for the fractional powers of nonlocal scalars we employ such as f J77H78|) . It should also be 
pointed out that the notion of energy for a nonlocal model is subtle, and more study of this 
issue is certainly required. If our class of models should prove to be unstable, it might be 
that the time scale is c/a ~ 6/H Q , which does not seem to pose a problem for galaxy and 
cluster dynamics. It might even be that the instability merely forces the weak fields back 
into the regime of general relativity which is stable. 



ACKNOWLEDGMENTS 

We have benefitted from discussions and correspondence with J. Bekenstein, S. McGaugh, 
M. Milgrom, R. Sanders, and C. Skordis. This work was partially supported by European 
Union grant INTERREG-IIIA, by NSF grant PHY-0855021, by the Institute for Fundamen- 
tal Theory at the University of Florida, as well as the ANR grant "THALES" . 



[1] C. M. Will, Living Rev. Rel. 9, 3 (2006), < ]http://www.livingreviews.org/lrr-2006-3[ >, 
[arXiv:gr-qc/0510072| ; K. Nakamura et al. [Particle Data Group], J. Phys. G 37 (2010) 075021. 

[2] G. Bertone (ed.), Particle dark matter: Observations, models and searches, (Cambridge Uni- 
versity Press, 2010); L. Bergstrom, New J. Phys. 11, 105006 (2009) [arXiv: 0903 .48491 [hep-ph]] . 

[3] M. Milgrom, Astrophys. J. 270, 365 (1983). 

[4] M. Milgrom, Astrophys. J. 270, 371 (1983). 



16 



M. Milgrom, Annals Phys. 229, 384 (1994) [arXiv:astro-ph/9303012~j . 

K. G. Begeman, A. H. Broeils and R. H. Sanders, Mon. Not. Roy. Astron. Soc. 249, 523 
(1991). 

E. F. Bell and R. S. de Jong, Astrophys. J. 550, 212 (2001) |arXiv:astro-ph / QQ 1 1493] . 

R. H. Sanders and S. S. McGaugh, Ann. Rev. Astron. Astrophys. 40, 263 (2002) 

[arXiv:astro-ph/0204521~] . 

R. B. Tully and J. R. Fisher, Astron. Astrophys. 54, 661 (1977). 

S. S. McGaugh, Phys. Rev. Lett. 106, 121303 (2011) |arXiv:1102.39l3l [astro-ph.CO]]. 

M. Kaplinghat and M. S. Turner, Astrophys. J. 569, L19 (2002) |arXiv:astro-ph/0107284~] . 

R. H. Sanders. larXiv:0806.2585l [astro-ph], 

R. H. Sanders, in The Dark Universe: Matter, Energy and Gravity, Proceedings of the Space 

Telescope Science Institute Symposium 15, ed. M. Livio (Cambridge University Press, 2003), 

p. 62, arXiv:astro-ph/0106558, 

R. H. Sanders. larXiv:1002.2765l [astro-ph.CO]. 

S. M. Faber and R. E. Jackson, Astrophys. J. 204, 668 (1976). 

S. S. McGaugh and W. J. G. de Blok, Astrophys. J. 499, 66 (1998) [arXiv:astro-ph/9801102] ; 
Astrophys. J. 499, 41 (1998) [arXiv:astro-ph/9801123l . 

W. J. G. de Blok and S. S. McGaugh, Astrophys. J. 508, 132 (1998) [arXiv:astro-ph/9805120] . 
X. Hernandez, M. A. Jimenez, and C. Allen, arXiv:1105.1873 [astro-ph. GA]. 
S. Dhital, A. A. West, K. G. Stassun and J. J. Bochanski, Astron. J. 139, 2566 (2010) 
|arXiv:1004.2755 [astro-ph.SR]]. 

A. Aguirre, J. Schaye and E. Quataert, Astrophys. J. 561, 550 (2001) 
|arXiv:astro-ph/0 105 184] . 

D. Clowe, A. Gonzalez and M. Markevitch, Astrophys. J. 604, 596 (2004) 
|arXiv:astro-ph/0312273~] . 

D. C. Dai, R. Matsuo and G. Starkman, Phys. Rev. D 78, 104004 (2008) |arXiv:0806.43l9l 
[gr-qc] ] • 

M. Feix, C. Fedeli and M. Bartelmann. larXiv:0707.0790l [astro-ph], 

G. W. Angus and S. S. McGaugh. larXiv:0704.038T1 [astro-ph]. 

J. D. Bekenstein, Phys. Rev. D 48, 3641 (1993) |arXiv:gr-qc/9211017| . 

J. D. Bekenstein and R. H. Sanders, Astrophys. J. 429, 480 (1994) |arXiv:astro-ph/9311062~] . 
R. H. Sanders, Astrophys. J. 480, 492 (1997) [arXiv:astro-ph/9612099] . 
J. D. Bekenstein, Phys. Rev. D 70, 083509 (2004) [Erratum-ibid. D 71, 069901 (2005)] 
[arXiv:astro-ph/0403694] . 

J. D. Bekenstein, PoS JHW2004, 012 (2005) |arXiv:astro-ph/0412652| . 



R. H. Sanders, Mon. Not. Roy. Astron. Soc. 363, 459 (2005) arXiv:astro-ph/0502222 1. 
C. Skordis, D. F. Mota, P. G. Ferreira and C. Boehm, Phys. Rev. Lett. 96, 011301 (2006) 
arXiv:astro-ph/0505519 j. 

C. Skordis, Phys. Rev. D 74, 103513 (2006) [arXiv:astro-ph/0511591]. 



S. Dodelson and M. Liguori, Phys. Rev. Lett. 97, 231301 (2006) |arXiv:astro-p h/0608602 1 . 
F. Bourliot, P. G. Ferreira, D. F. Mota and C. Skordis, Phys. Rev. D 75, 063508 (2007) 
|arXiv:astro-ph/0611255 j. 

T. G. Zlosnik, P. G. Ferreira and G. D. Starkman, Phys. Rev. D 77, 084010 (2008) 
|arXiv:071 1.0520 [astro-ph]]. 

C. Skordis, Class. Quant. Grav. 26, 143001 (2009) |arXiv:0903.3602l [astro-ph.CO]]. 
T. Clifton, P. G. Ferreira, A. Padilla and C. Skordis. larXiv: 1106. 24761 [astro-ph.CO], 



17 



[38] J.-P. Bruneton and G. Esposito-Farese, Phys. Rev. D 76, 124012 (2007) [ arXiv:0705. 4043 
[gr-qc]]. 

[39] I. Ferreras, N. E. Mavromatos, M. Sakellariadou, M. F. Yusaf, Phys. Rev. D80, 103506 (2009). 

|arXiv:0907. 14631 [astro-ph.GA]]. 
[40] N. E. Mavromatos, M. Sakellariadou, M. F. Yusaf, Phys. Rev. D79, 081301 (2009). 

|arXiv:0901 .39321 [astro-ph.GA]]. 
[41] M. A. Clayton (2001), |a?XTv:gr-qc/01041~03~} 

[42] M. D. Seifert, Phys. Rev. D 76, 064002 (2007) |arXiv:gr-qc/0703060| . 
[43] E. Sagi, Phys. Rev. D 80, 044032 (2009) [arXiv: 0905 .40011 [gr-ac]]. 

[44] J. D. Bekenstein, in Particle dark matter: Observations, models and searches, ed. G. Bertone, 

(Cambridge University Press, 2010). larXTv: 100 1.3876 1 [astro-ph.CO]. 
[45] T. Jacobson and D. Mattingly, Phys. Rev. D 64, 024028 (2001) [arXiv:gr-qc/0007031] . 
[46] T. Jacobson and D. Mattingly, Phys. Rev. D 70, 024003 (2004) |arXiv:gr-qc/0402005] . 
[47] C. Eling, T. Jacobson and D. Mattingly, |arXiv:gr-qc/0410001| 

[48] B. Z. Foster and T. Jacobson, Phys. Rev. D 73, 064015 (2006) | arXiv:gr-qc/0509083| . 
[49] T. Jacobson, PoS QG-PH, 020 (2007) [arXiv:0801.1547l [gr-ac]]. 

[50] A. Halle, Hong Sheng Zhao and Baojiu Liu, Astrophys. J. Suppl. 177, 1 (2008); Hong Sheng 

Zhao, Astrophys. J. 671, LI (2007). 
[51] E. Babichev, C. Deffayet and G. Esposito-Farese, Phys. Rev. D84, 061502 (2011). 

|arXiv:1106.2538l [gr-qc]]. 
[52] B. Li, H. Zhao, Phys. Rev. D80, 064007 (2009). [arXiv:0904.2935l [astro-ph.CO]]. 
[53] M. Milgrom, Phys. Rev. D 80, 123536 (2009) [arXiv:0912.0790l [gr-ac]]. 
[54] M. E. Soussa and R. P. Woodard, Phys. Lett. B 578, 253 (2004) [arXiv:astro-ph/0307358l . 
[55] R. Palais, Commun. Math. Phys. 69 (1979) 19. 

[56] C. G. Torre, AIP Conf. Proc. No. 1360 (AIP, New York, 2011), p. 63. [arXiv: 101 1.34291 [math- 

Ph]]- 

[57] N. C. Tsamis and R. P. Woodard. larXiv:1103.5134l [gr-ac], 

[58] N. C. Tsamis and R. P. Woodard, Phys. Rev. D 80, 083512 (2009) [ arXiv:0904.2368l [gr-qc]]; 
Phys. Rev. D 81, 103509 (2010) [arXiv: 1001 .49291 [gr-ac]]. 

[59] J. Schwinger, J. Math. Phys. 2, 407 (1961); K. T. Mahanthappa, Phys. Rev. 126, 329 (1962); 
P. M. Bakshi and K. T. Mahanthappa, J. Math. Phys. 4, 1 (1963); J. Math. Phys. 4, 12 
(1963); L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1965); K. C. Chou, Z. B. Su, B. L. Hao 
and L. Yu, Phys. Rept. 118, 1 (1985); R. D. Jordan, Phys. Rev. D 33, 444 (1986); E. Calzetta 
and B. L. Hu, Phys. Rev. D 35, 495 (1987). 

[60] B. S. DeWitt and R. W. Brehme, Annals Phys. 9, 220 (1960). 

[61] S. Deser, Phys. Lett. B 479, 315 (2000) [ arXiv:hep-th/9911129| ; Nucl. Phys. Proc. Suppl. 88, 
204 (2000). 

[62] S. Park and R. P. Woodard, Gen. Rel. Grav. 42, 2765 (2010) |arXiv:0910.4756l [gr-ac]]. 
[63] N. C. Tsamis and R. P. Woodard, Annals Phys. 267, 145 (1998) [arXiv:hep-ph/9712331] . 
[64] M. E. Soussa and R. P. Woodard, Class. Quant. Grav. 20, 2737 (2003) 

|arXiv:astro-ph/0302030] . 
[65] G. L. Page, D. S. Dixon and J. F. Wallin, Astrophys. J. 642, 606 (2006); J. F. Wallin, 

D. S. Dixon and G. L. Page, Astrophys. J. 666, 1296 (2007) [arXivj0705j4Q8] [astro-ph]]. 
[66] X. Wu, H. Zhao, B. Famaey, G. Gentile, O. Tiret, F. Combes, G. W. Angus, A. Robin, 

Astrophys. J. 665, L101 (2007) [arXiv:0706.3703l [astro-ph]]. 



18 



